Splitting of the Raman 2D band of graphene subjected to strain 
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The Raman 2D-band -important for the analysis of graphene- shows a splitting for uniaxial strain. The 
splitting depends on the strength and direction of the applied strain and on the polarization of the incident 
and outgoing light. We expand the double-resonance Raman model in order to explain the strain direction 
dependence and the polarization dependence of the splitting. The analysis of this splitting gives new insight into 
the origin of the 2D-band. Our prediction of the strain direction and polarization dependence agrees well with 
recent experiments. 
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The discovery of graphene in 2004 has led to strong re- 
search activities in the last years QJJ] . Graphene has been shown 
to possess unique material properties. In graphene the quan- 
tum Hall effect was observed at room temperature QS]. Be- 
cause of the specific band structure near the Fermi level, the 
electrons behave like massless Dirac fermions and mimic rel- 
ativistic particles with zero rest mass and with an effective 
'speed of light' c' 10 ■ 10 6 m/s QH. The high electron 
mobility makes graphene a promising candidate to serve as 
building block for micro-electronics. 

Detecting and characterizing graphene still is an exper- 
imentally difficult task, in particular, finding its crystallo- 
graphic orientation. One method to specify graphene layers 
is Raman spectroscopy. This non-destructive method can be 
used at ambient conditions. It allows to distinguish between 
single-layer, bilayer, and multi-layers of graphene (5]. The 
quality of the graphene flakes can be estimated from Raman 
features jfj. Recently, uniaxially strained graphene has been 
investigated by Raman spectroscopy JJS]. The strain softens 
the frequency of the optical phonon branches. Uniaxial strain 
also reduces the hexagonal symmetry of the system. This 
symmetry breaking leads to new effects, e.g. , the orientation 
of graphene can be determined from the polarization depen- 
dence of the split G-bandj8l|9t]. Because of the doping depen- 
dence of the G-band-frequency fioll and the higher shift rate 
of the 2_D-mode the latter is more commonly used to quantify 
strain in graphene. The shift of the 2D-mode is currently un- 
der discussion and usually only a single value for all measure- 
ments is given. S EH] However, recent experiments show 
that the shift of the 2£>-mode in addition depends on the crys- 
tallographic orientation and the polarization of the incident 
and outgoing light during the Raman experiments . 1 1 21 1 1 311 In 
addition, a splitting of the 2£>-mode has been observed. To 
understand these new results a comprehensive theoretical de- 
scription is needed to quantify the shift of the 2£>-mode under 
strain and that takes into account the polarization and the crys- 
tallographic orientation. 

Here we expand the established model of the double- 
resonance in graphene. We explicitly calculate the resonant 
Raman cross section which is used to identify the contribut- 
ing phonons and to estimate a polarization dependence. The 
effect of small strain on the peak positions of the 2D mode of 



graphene is investigated on the level of density-functional the- 
ory. The strain is applied along arbitrary directions. We are 
able to explain the splitting, the shift rates, the polarization 
effects and the crystallographic orientation dependence. The 
observed effects do not solely come from a movement of the 
Dirac cones or their deformation, but rather from an orienta- 
tion dependent softening of the involved phonon TO branch. 
The large magnitude of the splitting can only be explained 
with phonons from the branches between the T and A' -points 
that contribute to the double-resonant 2£>-mode. 
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FIG. 1. Double-resonance mechanism (left: outer process, right: in- 
ner process) plotted into the band structure scheme of graphene along 
KMK' . Below is the mechanism plotted in the Brillouin zone. 

In Fig.Q]we show the scheme of the double-resonance prin- 
ciple for the 2Z?-mode[ 14] . An electron-hole pair is created by 
an incoming photon with energy Ei. The electron or hole is 
scattered inelastically by a phonon with a wavevector qi ^ 
into a real state. From there it scattered back inelastically by 
another phonon q2 = —qi, where it recombines with an en- 
ergy E = Ei — 2hu(qi). The difference to the 2D-mode is 
that for the _D-mode, one inelastic scattering process is re- 
placed by an elastic scattering at a defect. It is generally 
accepted, that the contributing phonons stem from the fully- 
symmetric TO-derived branch] 15]. 

In the following we first use analytical expressions of the 
electronic and vibrational band structure to calculate Raman 
spectra for different scattering configurations. With these re- 
sults we identify the dominant processes and motivate their 
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extraction from density functional theory (DFT) calculations 
under strain. The cross section |AT2/,io| 2 of a resonant Raman 



process can be calculated bv lU6H 
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This is the formula for the D-mode. The first term denotes 
incoming resonance, and the second outgoing resonance. E xy 
is the energy difference between the electronic states x and y. 
Ei and Ei — hw are the energies of the incoming and outgo- 
ing photon, respectively. M xy is the matrix element for the 
scattering over the intermediate states x and y, and 7 is the 
broadening parameter of the electronic transition. 

To investigate which processes contribute to the Z?-mode 
we evaluate Eq. £[). The sum is converted into an integral and 
the matrix elements are taken constant, although they may 
depend on the electronic wavevector. The electronic bands 
are from a tight-binding approach fl7ll . where we use the val- 
ues 70 = 2.8 eV and 70 ' = OeV. The TO phonon dispersion 
around the K -point is taken from the analytical expression de- 
rived from recent inelastic X-ray measurements [18]. 

In a first approximation we integrate along the line that 
connects two neighboring K-points [see Fig.|2a)]. We have 
to take special care in this one-dimensional integration as it 
leads to erroneous results: in the one-dimensional integrations 
phonons with q w K do not interfere destructively and give a 
major contribution (3. Still, the one-dimensional integration 
gives a good idea for identifying the contributing processes. 
We find contributions from phonons from the K — M and 
r — K directions. The Raman matrix element |if2/,io| as a 
function of phonon frequency is plotted in Fig.[2|b). Here, a 
small broadening parameter of 7 = 0.005 eV was used. The 
two peaks of the blue and green plot correspond to incoming 
and outgoing resonance. 

In a second approximation we now look at the integration 
in two dimensions. The integration areas are highlighted as 
shaded areas in Fig. |2ja), where scattering into all neighbor- 
ing AT-valleys is considered. The integral is evaluated with 
the Monte-Carlo method provided in the Mathematica pack- 
age (the method ' AdaptiveMonteCarlo' was used with at least 
5 x 10 6 points) fl^l . Now, the AT -phonon contribution (not 
shown) is cancelled out due to destructive interference, in 
agreement with Ref. |fl5ll . Comparing the result in Fig. |2|b) 
with the one-dimensional integration we find similar features. 
We thus can conclude that phonons from the K — AI direc- 
tion and phonons from the T — K direction do contribute to 
the D-mode in graphene. We will refer to them as outer pro- 
cesses (phonons from K — M, see Fig.[T|i and inner processes 
(phonons from T — K). Using a larger value of 7 = 0.05 eV, 
the distinct features merge and a single broad D-band arises 
[dotted curve in Fig.|2b)]. 

In a perfect graphene flake the D-band does not depend on 
the polarization of the incident light. The polarization of the 
incoming laser affects the absorption in fc-space. Electrons 
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FIG. 2. (a) The two-dimensional integration area plotted in the elec- 
tronic energy contours. The electron absorption happens near "ai" 
and the electron is scattered to regions near "bi". (b) Raman matrix 
element | A*2 /.10 1 f° r several different parameter sets. Integration is 
either performed along lines connecting neighboring A'-points (de- 
noted by '1-dim.') or two-dimensional within the highlighted area 
('2-dim.'). The similar features in both types of integration verifies 
that contributions from the inner processes are visible in the full cal- 
culation (two-dimensional integration) as well. At higher values of 
7, single contributions are not distinguishable and a single D-peak is 
visible. 



with a wavevector parallel to the polarization of a photon can 
not be excited during an absorption process. We include the 
effect of polarization by using the \p x fc| 2 -dependence of the 
absorption, where p denotes the polarization of the incident 
light and k denotes the wavevector of the excited electron ic . 
This leads to areas in the Brillouin zone that give zero con- 
tribution to the Raman cross section. This is illustrated in 
Fig. 0a) for "vertical" polarization indicated by the arrow p. 
As expected and known from experiments, polarization of the 
incoming light has no influence on the D-band for unstrained 
graphene in Fig. 0b). Now we simulate tensile strain, by com- 
pressing the band structure in strain-direction. This is indi- 
cated in Fig. 0a) for strain in the zigzag direction. A strain 
of 9% was used. This is just a proof of concept, as we have 
not included the phonon softening, and used the unstrained 
phonon dispersion. The resulting Raman cross sections are 
shown in Fig. 0b). The labels parallel and perpendicular are 
w.r.t. the strain direction. Depending on the polarization, the 
£>-band has a different energy shift. 

Using DFT calculations of the phonon dispersion and the 
electronic band structure we now make quantitative predic- 
tions on the peak positions of the 2D band for several scatter- 
ing configurations. We have calculated the electronic bands 
and phonon dispersion of graphene under strain. These 
calculations were done with the linear-response DFT code 
QUANTUM-ESPRESSOIMlISI- Uniaxial strain was applied 
via the two-dimensional strain tensor e = ((e, 0), (0, — e^)), 
where v is the Poisson's ratio. Strain in arbitrary directions 
has been applied through rotating the strain tensor j23ll . The 
strain (tensile and compressive) is applied in several directions 
between the zigzag and armchair direction for strain values up 
to 3 % J22ll23ll However, the largest differences are found be- 
tween zigzag and armchair direction. With this data we can 
now perform a mapping of the expected 2£>-band frequencies, 
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FIG. 3. (a) Modelling of strain by compressing the band structure. 
To account for the anisotropic absorption depending on the polar- 
ization, parts of the integration area have been omitted, (b) Raman 
matrix element |iv2/,io| for strained and unstrained graphene with 
polarization parallel and perpendicular to the strain direction. Only 
when strain is applied, a polarization dependent shift is observed. 



using the scattering processes shown in Fig.Q] The scattering 
of the electrons is considered to be elastic for identifying the 
contributing phonon's wavevector (this corresponds to hori- 
zontal arrows in Fig . [TJ . 

Depending on which direction in reciprocal space the elec- 
tron is scattered, different resonance conditions apply. This 
is a consequence of the symmetry-breaking induced by the 
strain. J3 In addition, as mentioned above, when the light 
is polarized, some areas in the reciprocal space do not have 
excited electrons that contribute to the 2D-band. With this 
argument the most dominant processes are those that are per- 
pendicular to the polarization of the incoming light. Fig. [4] 
now show the obtained 2D frequencies for inner and outer 
processes (as depicted in Fig. [TJ, for strain in armchair (AC) 
and zigzag (ZZ) direction, and for two different laser ener- 
gies (1.5 eV and 2.4 eV). The corresponding processes are 
illustrated in Fig. HJa). The notation K\ means along K\, 
the AT-points are numbered clockwise, beginning with the top 
right as number 1. As can be seen, the splitting is larger for 
ZZ strain. Also the polarization dependent shift is larger for 
ZZ strain. Both, inner and outer processes lead to a polariza- 
tion dependent splitting. However, the inner processes lead 
to a larger splitting. In Fig. |2b) we show the phonon disper- 
sion along two different paths in the reciprocal space. The TO 
branch between Y and K shows a large splitting depending 
on the chosen path. The softening is stronger in the direction 
perpendicular to the strain direction. This can be understood 
looking at the displacement patterns of this mode: the dis- 
placement is parallel to the strain direction and the involved 
atomic bonds are weakened due to the strain. The same ar- 
gument is responsible for the stronger softening of the G~- 
mode at the zone-centerJUIg, 24]. Between K and M the dis- 
placement patterns of the TO branch are not parallel anymore 
and have lost their TO character and thus are not affected so 
strongly by the strain direction. 

Biaxial strain was also applied to the system. In this case, 
the symmetry is conserved, there is no splitting. For the G- 
band we obtain a shift rate of 60 cm -1 / %. This is in excellent 
agreement with the value of 63 cm -1 / % obtained in Ref. M25I1 



. For the 2D-band we find a shift rate of 135 cm~V% for both 
inner and outer processes. 

With the above method of frequency mapping we cannot 
determine the details of intensity and lineshape of the D and 
2£>-bands. Simple arguments like looking at the electronic 
density of states do not take into account the phonon density 
of states. Sometimes it is argued that the outer processes 
only contribute to the 2£>-mode because of the trigonal warp- 
ing effect and a higher density of states 1261 271: The trigonal 
warping effect becomes stronger for larger excitation energies, 
however for optical transition energies <2 eV the equi-energy 
contours are still more or less round. In addition, the phonon 
dispersion of the TO-derived mode near the A'-point was a 
long time believed to show a strong trigonal warping. This 
was recently resolved by Griineis et al. performing inelastic x- 
ray experiments revealing almost-round equi-energy phonons 
of the TO-branch near the A'-point[ 18] with almost no trigo- 
nal warping. 
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FIG. 4. Mapping of the Raman 2D contributions for inner and outer 
processes for the laser energies 1.5 eV (a) and 2.4 eV (b). Circles 
(triangles) correspond to strain in zigzag (armchair) directions. Open 
(closed) symbols denote parallel (perpendicular) polarization w.r.t. 
the strain direction. The data points are connected with lines that 
denote inner (dotted) and outer (dashed) processes. 

In summary, we have presented an in-depth analysis of 
the evolution of the D and 2£>-bands in uniaxially strained 
graphene. Depending on polarization and strain direction, 
different shifts of the Raman-active D and 2£>-band are pre- 
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FIG. 5. (a) Dominant processes for perpendicular (solid lines) and 
parallel polarizations (dashed lines) for strain in ZZ and AC di- 
rection, (b) Optical phonon branch along TKM in graphene for 
unstrained (blue, solid line) and 2% strain in ZZ direction (green, 
dashed-dotted line and red, dashed line). The corresponding paths in 
the Brillouin zone are shown in the inset for an exaggerated strain 
of 30%. The dark grey (blue) arrow indicates the outer process for 
a laser energy of 2.4 eV, the grey (green) arrow indicates the inner 
process. Note that in absolute units of A -1 the three branches are 
slightly different due to strain (see inset). In the main figure we 
scaled the wavevector axis such that all K and M points are at the 
same position. 



dieted. Strain in zigzag directions leads to a larger splitting 
than strain in armchair directions. Our predictions are in 
excellen t ag reement with experiments on strained graphene 
sample s i 121 1 1311 . As the strain-induced shift is used to deter- 
mine strain in graphene our results are relevant for the inter- 
pretation of experimental data. 
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